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ABSTRACT 

To use galaxy clusters as a cosmological probe, it is important to account for 
their triaxiality. Assuming that the triaxial shapes of galaxy clusters are induced 
by the tidal interaction with the surrounding matter, Lee and Kang recently 
developed a reconstruction algorithm for the measurement of the axial ratio of 
a triaxial cluster. We examine the validity of this reconstruction algorithm by 
performing an observational test of it with the Virgo cluster as a target. We 
first modify the LK06 algorithm by incorporating the two dimensional projec- 
tion effect. Then, we analyze the 1275 member galaxies from the Virgo Cluster 
Catalogue and find the projected direction of the Virgo cluster major axis by 
measuring the anisotropy in the spatial distribution of the member galaxies in 
the two dimensional projected plane. Applying the modified reconstruction al- 
gorithm to the analyzed data, we find that the axial ratio of the triaxial Virgo 
cluster is (1: 0.54 : 0.73). This result is consistent with the recent observational 
report from the Virgo Cluster Survey, proving the robustness of the reconstruc- 
tion algorithm. It is also found that at the inner radii the shape tends to be more 
like prolate. We discuss the possible effect of the Virgo cluster triaxiality on the 
mass estimation. 

Subject headings: cosmology: theory — large-scale structure of universe 

1. INTRODUCTION 

Galaxy clusters provide one of the most powerful tools to constrain the key cosmological 
parameters. In the era of precision cosmology, it is important to determine their mass as 
accurately as possible before using them as a cosmological probe. Any kind of simplified 
assumption about the properties of galaxy clusters could cause substantial systematics in 
the mass estimation. The triaxial shapes of galaxy clusters are one of such properties. 
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It has been long known both obseryationaUy and numerically that the galaxy clusters 



are n oticeably triaxial (e.g., iFrenk et al.lll988l ; IWestlll989l ; iPlionis et al.lll99ll ; IWarren et al. 



1992). While plenty of efforts have been already made to take into account 


;he triaxial s 


lapes 


of galaxy clusters (Jing & Suto! 20021; 


Fox & Pen 2001: Suwa et al. 20031: 


Lee & Suto! 


2004: 


Kasun & Evrard 20051; Hopkins et al.l 


20051: Lee et al. 20051: Smith & Watts 


12005c Havashi et al. 


2007h . previous studies have been largely focused on the statistical treatment of cluster tri- 



axiality. For the measurement of the gas mass fraction of galaxy clusters that can pro 
vide powerful cons t raints on the den s ity param e ter and the dark energy equation of state 



rtwhite et al.lll993l: iLubin et al.lll996l : [Cenlll997l : [Evrardlll997l : ICooravi nmS . lLaroque et al. 



2006 



Ferramacho fc Blanchardl 120071 ) . however, it is necessary to deal with the individual 



clusters and their triaxial shapes. 



The standard picture based on the cosmic web theory (IBardeen et al.lll986l ; iBond et al. 



19961 ) explains that galaxy clusters are rare events corresponding to the local maxima of 
the initial density field and form at the dense nodes of the local filaments in the cosmic 
web through tidal interactions with the surrounding matter. The tidal effect from the sur- 
rounding matter results in the deviation of cluster shapes from spherical symmetry as well 



as the pr eferentia 



alignments of cluster galaxies (or halos) with the major axes of their host 



clusters ( 


3ine;eeli 


1982; 


Strub 


e k, Peebles 


1985; 


Hop 


tins et al. 


2005; 


Kasun & Evrard 


2005 


Lee et al. 


2005: 


Aleood et al. 


2006; 


Altav et al. 


2006; 


Paz et al 


2006 


)■ 



In the frame of this standard scenario, iLee &; Kang I (120061 . hereafter, LK06) have re- 
cently developed an analytic algorithm to reconstruct the triaxial shapes of individual clus- 
ters. The key concept of the LK06 algorithm is that the two axial ratios of a triaxial cluster 
are related to the eigenvalues of the local tidal shear tensor. By measuring the spatial align- 
ment of the cluster galaxies with the major axis of their host cluster, one can determine the 
eigenvalues of the local tidal tensor, which will in turn yields the two axial ratios of a triaxial 
cluster. 

Testing their analytic model against high-resolution N-body simulations, LK06 have 
shown that their algorithm works well within 20% errors. Now that the LK06 algorithm is 
known to work in principle, it is time to test the algorithm against observations. Our goal 
here is to apply the LK06 algorithm to real observational data and examine its validity in 
practice. Here, we use the Virgo cluster as a tar get, whose triaxial shape has been very 
recently measured observationally (jMei et al.l 120071 ). 



The organization of this paper is as follows. In §2, a brief overview of the LK06 algorithm 
is provided and how to incorporate the two dimensional projection effect into the algorithm is 
explained. In §3, the Virgo cluster data are analyzed and its triaxial shapes are reconstructed 
using the LK06 algorithm. In §4, the results are summarized, and the advantages and the 
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caveats of our model are discussed. 



2. THEORETICAL MODEL 

2.1. Overview of the LK06 Algorithm 

According to the LK06 algorithm, the cluster triaxial shape originates from its tidal 
interaction with the surrounding matter distribution. This assumption has been verified from 
high-resolution N-body simulation which demonstrated clearly that the tidal field elongates 



the cluster shapes (e.g.. lAltay et al.ll2006l . and references therein). LK06 has shown that the 
two axial ratios of a triaxial cluster are related to the three eigenvalues of the local tidal 
tensor defined as the second derivative of the gravitational potential: 



l / 1 \ \ 1/2 /-, \ \ 1/2 

b I 1 — A 2 \ a i 1 — Ai x 



c \ 1 — Aq / c V 1 — A 



where {a, b, c} (with a < b < c)are the three principal axis lengths of a cluster and {Ai, A 2 , A 3 } 
(with Ai > A 2 > A3) are the three eigenvalues of the local tidal tensor, T. According to this 
formula, one can estimate the cluster axial ratios if Ai, A2, A3 is fixed by fitting the probability 
density distribution, ^(cos^) analytically to the observational data. LK06 suggested that 
the preferential locations of the cluster galaxies near the cluster major axes given the local 
tidal tensor be described as 

1 — s - - 

(xiXj\T) = — — 8ij + sT ik T kj . (2) 

where x = is the unit position vector of a cluster galaxy, T = T/|T| is the unit tidal 
shear tensor, and s G [—1,1] is the correlation parameter that represents the correlation 
strength between x and T. Under the assumption that the minor axis of the tidal shear 
tensor is in the direction of the cluster major axis, equation ([2]) basically describes the 
alignment between the position of a cluster galaxy and the major axis of its host cluster. If 
s = — 1, there is a maximum alignment. If s = 1, there is a maximum anti-alignment. The 
case of s = corresponds to no alignment. 

Let 9 3 d be the angle between the host cluster major axis and the galaxy position vector. 
Under the assumption that the cluster major axis is in the direction of the minor principal 
axis of the tidal shear tensor, The probability density distribution of cos 9 3 d was derived by 
LK6 as 

p{cos6 3d ) = . M~ x ■ %)-!, (3) 
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where is an azimuthal angle of x measured in an arbitrary coordinate system. Here, 
covariance matrix M is defined as M = (xiXj\T). Note that equation ([3]) holds good for any 
arbitrary coordinate system in which the tidal shear tensor is not necessarily diagonal. 

In the principal axis frame of the tidal tensor, equation (j3J) can be expressed only in 
terms of the three eigenvalues of the tidal shear tensor as 

p(cos 6) = J_f[(l- S + 3 S A7)- 1/2 



2tt 

i=i 



-3/2 



/ sin 6 cos d> sin 2 9 sin 2 6 cos 2 6 . 

£- + £- + — dcp, 4 

1 - s + 3sA 2 l-s + 3s\ 2 2 l-s + 3s\ 2 3 ' 

where {Aj}f =1 are the eigenvalues of T, related to {Aj}f =1 as 

Ai = s 6X s , A 2 = s 4 3 2 s , A 3 = s s , (5) 

Ai + A2 + A3 Ai + A2 + A3 Ai + A2 + A3 



where 5 C ~ 1.68 is the linear density threshold for a dark halo (jEke et al.lll996l ) satisfying 
the following constraint of 5 C = Yli=i ^i- 

The key concept of LK06 algorithm is that by fitting equation (j3j) to the observed 
probability distribution, one can find the best-fit values of Ai, A2 and s, and then determine 
the cluster axial ratios using equation (TjQ). Although LK06 algorithm allows us in principle 
to reconstruct the three dimensional intrinsic triaxial shapes of individual clusters, it is 
restricted to the cases where the informations on the three dimensional positions of the 
cluster galaxies in the cluster principal axis frame are given. Unfortunately, for most clusters, 
these informations are not available but only two dimensional projected images of clusters. 

In §2.2, we attempt to modify the LK06 algorithm in order to incorporate the two 
dimensional projection effect. 



2.2. Projection Effect 

Let us suppose that the position vectors of the cluster galaxies are all projected along 
the line of sight direction onto the plane of a sky. Unless the major axis of the host cluster is 
perfectly aligned with the line-of-sight direction to the cluster center, one would expect that 
the projected position vectors of the cluster galaxies should show a tendency to be aligned 
with the projected major principal axes. 

Let 02d be the angle between the projected cluster major axis and galaxy position vector 
in the plane of sky. The probability distribution can be calculated by integrating equation 
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j3]) along the line of sight as 

p(cos 9 2d ) 



1 Z7T 



(6) 



where (£3) is now chosen to be in the direction of the line of sight. In other words, we 
consider a certain Cartesian coordinate system in which the x 3 direction is parallel to the 
line of the sight to the cluster center of mass. Note that in this coordinate system the tidal 
tensor is not necessarily diagonal. 



Through the similarity transformation 

f = R* A T R, 



(7) 



with 
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one can express the nondiagonal u nit tidal tens or T in terms of its eigenvalues. Here, the 
rotation matrix R has the form of (jBinneylll985l ) 



- sin ip — cos ip cos £ cos ip sin £ 
R = I cos ip — sin ip cos £ sin ip sin £ 
sin £ cos £ 



(9) 



where {£,ip) is the polar coordinate of the line-of-sight direction in the principal axis frame 
of the cluster. 

Through equations ©-Q, we finally final an analytic expression for the probability 
density distribution of cos92d in terms of {Aj}| =1 and the correlation parameter s: 



p(cos 9 2 d) 
with the factor Q defined as 



:i + *) 2 (i-2s) 

2tt 



Q- 3 / 2 d 



X3, 



(10) 



-1 



Q 



[(1 + s) 2 - 3s{l + s)]{A x Xi + A 2 \ 2 + A 3 A 3 )(1 



X 



3s(l + s)](C 1 \ 1 + C 2 \ 2 + C 3 \ 3 )x 



+ [(1 



+6s(l + s)(B 1 X 1 + B 2 X 2 + B 3 X 3 )x 3 Jl 



x 
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where the coefficients {Ai}^ =1 , {i? i }? =1 , {Cj}? =1 are given as 

A\ = cos ip cos 2 d( cos ip cos &2d — 2 sin ip cos £ sin 9 2d) + sin 2 ^(sin 2 ip + sin 2 £ cos 2 ip) , 



- 6- 



A 2 = sin Tp cos 9 2 d (sin ip cos 9 2 d + 2 cos ^ cos £ sin 6*2d) + sin 2 9 2 d (cos 2 ?/> + sin 2 £ sin 2 -0) , 

A 3 = cos 2 6> 2(J + cos 2 £ sin 2 # 2 d, 

£>i = — sin £ cos ip (cos £ sin cos ip 4~ sin cos 9 2( i) j 

£> 2 = — sin£ sin-0(cos£sin ^sin^ — cosip cos 9 2d ), 

B 3 = cos £ sin ^ sin 9 2 d, 

C\ = sin 2 + cos 2 ^ cos 2 -0, 

C 2 = cos 2 ip + cos 2 £ sin 2 ■0, 

C 3 = sin 2 e- 

With this new modified algorithm in the two dimensional space, we can reconstruct the 
intrinsic shape of a triaxial cluster halo from the observed two dimensional image. 



APPLICATION TO THE VIRGO CLUSTER HALO 



3.1. Observational Data and Analysis 

The Virgo cluster is the nearest ric hly populated cluster of galaxies whose prope rties 
has been studied fruitfully for long (e.g., Bohringer et al. 1994 ; West fc Blakeslee I2OOO, and 



refere nces therein). It is known to have approximately 1275 member galaxies (IBinggeli et al. 



19851 ) and located at a distance of approximately 16.1 Mpc from us (ITonry et al.l 120011 ) with 
the major axis inclined at an angle of approximately 10° with respect to the line of sight 
rtWest fc Blakeslee! l2000h . 



Near the center of the Virgo cluster is located the large ellipticity galaxy M87 (or Virgo 
A). The major axis of the Virgo cluste r is found to be inclined a pproximately 10° with respect 
to the line of sight direction to M87 (jWest fe Blakesled 120001 ). We use data from the Virgo 
Cluster Catalogue (IBinggeli et al.lll985l ) which compiles the equatorial coordinates of total 
1275 member galaxies. 



3.2. Coordinate Transformation and Projection Effect 

Let r and (a, 5) represent the three dimensional distance to a member galaxy and its 
equatorial coordinates, respectively. Under the assumption that the position of M87 is the 
center of mass of the Virgo cluster, the Cartesian coordinate of a member galaxy in the 
center of mass frame can be written as 

xi = r cos a cos 5 — Xi va 
x 2 = r sin a cos 5 — x 2va 
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x 3 = rsm5-x 3v& (12) 

where (xi vs ,,x 2va ,,x 3va ) represents the position of M87. The equatorial coordinates of M87 
is measured to be a va = 187.71° and <5 va = 12.39°. The dista nce to Virgo A, r m , from us is 



also known to be approximately 16.1Mpc (ITonry et al.ll200ll ; SBF survey). Thus,we have a 
full information on (xi va , X2 va ,,x 3va ). 

Now let us consider a coordinate system where the third axis is in in the direction of the 
line of sight to M87. Let (u, <p) be the spherical polar coordinate of the member galaxy 
in this coordinate system. It can be found through coordinate transformation as 

va. 

COS -0 cos £ va cos ip — simjj 

va. 

'Q^ va sin^ va COS £ va sin Ipya COS^va ) ( %2 ) (13) 

cos£ 

va 

— sin£ 

va 









H 







-sinVVa ^ 






COS V'va 


II 


^2 


o ) 






™ = ir/2 - 


^va- 



where (£ va , ^a) is the polar coordinate of M87, so that £, 

In the plane of sky projected along the line of sight direction, the spherical polar coordi- 
nates of the member galaxy can be regarded as the Cartesian coordinates as r — >• 0, $ — > xi, 
and (j) — > yi- Basically, it represents a two dimensional projected position of a Virgo cluster 
member galaxy in the plane of sky with the position of M87 as a center. 

For each member galaxy from the Virgo Cluster Catalogue, we have determined the two 
dimensional projected position (xi,yi) using the given equatorial coordinates. 



3.3. Virgo Cluster Reconstruction 

To measure the alignments between the positions of the Virgo cluster galaxies and the 
projected major axis of the Virgo cluster and compare the distribution of the alignment 
angles with the analytic model (eq. [6]), it is necessary to find the direction of the projected 
major axis in the coordinate system of (xu,X2i)- Equivalently, it is necessary to have an 
information on the polar coordinates of the line of sight, (£, ip) with respect to the three 
dimensional principal axis of the Virgo cluster. 



The seminal paper of IWest fc Blakesled (120001 ) provides us with the information on £ 



i.e., the angle between the three dimensional major axis of the Virgo cluster and the line of 
sight direction as approximately 10°. However, we still need the azimuthal angle ip of the 
major axis of the Virgo cluster. 

To find the azimuthal angle ip, we first let £ = 10° and ip = in the analytic model 
(eq.[6]). It implies that we choose a certain Cartesian coordinates (xi p ,X2 P ) where the angle 
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ip vanishes. Then, we transform the coordinate system of (xu,X2i) into this new coordinate 
system as 

( Xlp \ = ( cos ^ sin ^ \( Xu \ (U) 
V X 2 P J V -sin^ cosip J \x 2 i J ' 

Here, note that (cos if), simp) corresponds to the projected major axis of the Virgo cluster in 
the {xu,X2i) coordinate system. It is expected that in this new (x\ p , X2 P )-coordinate system 
the the observationally measured distribution should fit the analytic model (eq.[6] best. 

For a given if>, we measure the alignment angle between the projected major axis and 
position vector of each Virgo cluster galaxy as 

cos#2ii = %i P cosip + X2 P sm if), (15) 

where x p = x p /|x p |. Then, by counting the number of galaxies galaxy's number density 
as a function of cos 9d, we can derive the probability distribution, p(cos9d). We fit this 
observational distribution with the analytic model, adjusting the values of Ai, A2 and s 
through ^-minimization. We repeart the whole process varying the value of if), and seek 
for the value of if) which yields the smallest x 2 value. As a final step, we determine the 
corresponding best-fit values of Ai and A2. 

Finally, we find the axial ratios of the Virgo cluster to be a/c = 0.54 and b/c = 0.73 
by using equation ([T]) with the constraint of 5 C = Y^=i ^t- The best-fit value of s is also 
determined to be —0.25, indicating that the positions of the cluster galaxies are indeed 
aligned with the major axes of the Virgo cluster. Figure [1] plots the two dimensional 
projected positions of the Virgo cluster galaxies in the (xi p , X2 P )-coordinate system. The 
arrow represents the direction of the projected major axis determined from the chosen value 

of if). 

Figure [2] plots the probability distribution of the alignment angles for four different cases 
of if). In each panel, the histogram with Poisson errors is the observational distribution while 
the solid line represents the analytic fitting model. The dotted line stands for the case of no 
alignment at all. The top left panel corresponds to the finally chosen value of if> for which 
the analytic fitting model and the observational result agree with each other best. The other 
three panels show the three exemplary cases of if) for which the analytic fitting model and 
the observational result do not agree with each other well with relatively high \ 2 value. 

To investigate whether the reconstructed axial ratios of the Virgo cluster changes with 
the radius from the center, we introduce a cut-off radius, -R cu t, and remeasure the axial ratios 
using only those galaxies located within R cnt from the center in the (x\ p , X2 P )- corrodinate 
system. We repeat the same process but using four different values of i? cu t- Table [T] lists the 
values of the resulting best-fit axial ratios a/c and b/c, and the best-fit correlation parameter 
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for the four different cut-off radii of R cu t- As can be seen, at the inner radius smaller than 
the maximu m one, the shape ten ds to be more prolate-like, consistent with recent numerical 
report (e.g.. lHayashi et al.ll2007l ). Note also that the value of s is consistently —0.25, which 
implies the strength of the tidal interaction with surrounding matter is consistent. 

Figure [3] plots the probability distributions for the four different cases of -R cut . The top 
left panel shows the case of maximum cut-off radius. As can be seen, the agreements between 
the analytic fitting model and the observational result are quite good even at inner radii. 



4. SUMMARY AND DISCUSSION 

We have modified the cluster reconstruction algorithm which was originally developed 



by lLee fc Kang I (120061 ) to apply it to the two dimensional projected images of galaxy clusters 
in practice. Assuming that unless the cluster major axes are in the line-of-sight direction, 
we have found that the alignments between the galaxy positions and the projected major 
axes can be used to reconstruct the two axial ratios of the triaxial clusters. We have applied 
the modified algorithm to the observational data of the Virgo cluster and shown that the 
reconstructed axial ratios of (1 : 0.54 : 0.73) a re in good agreement with the recent report 



from the Virgo Cluster survey (IMei et al.ll2007l ). which proves the validity and usefulness of 
our method in practice. 

Now that the Virgo cluster is found to be triaxial, let us discuss on the triaxiality 
effect on the mass estimation. For simplification, let us assume that the Virgo cluster has a 
uniform density. Then, the mass of the triaxial Virgo cluster with the axial ratios given as 
(1 : 0.54 : 0.73) is estimated to be 2.6 times larger than the spherical case since the spherical 
radius is close to the minor axis length of the Virgo cluster since the major axis of the Virgo 
cluster is very well aligned with the line of sight direction. Therefore, it can cause maximum 
~ 50% errors to neglect the Virgo cluster triaxiality. 

The most prominent merit of our method over the previous one is that it reconstructs 
directly the intrinsic three dimensional structures of the underlying triaxial dark matter halo 
using the fact that the cluster triaxiality originated from the tidal interaction. Convention- 
ally, the triaxial shape of a cluster is found through calculating its inertia momentum tensor. 
This conventional method, however, is unlikely to yield the intrinsic shape of the underlying 
dark matter halo unless the target cluster is a well relaxed system. In addition, our method 
does not resort to any simplified assumption like the axis-symmetry or the alignment with 
the line-of-sight and etc. 



Yet, it is worth mentioning here a couple of limitations of our method. First, it assumes 
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that the major axes of the clusters are not aligned with the line-of-sight direction so that 
the alignments between the galaxy positions and the projected major axes can be measured. 

Second, for the Virgo cluster, the crucial information on the angle £ between the thre e 



dimensional major axis and the line of sight has been already given IWest fc Blakesled (120001 ) . 
Which simplifies the whole of our method since it was only if) that has to be determined. 
For most clusters, however, this information is not given, so that both the values of £ and ip 
have to be determined through fitting before finding the axial ratios. 

Third, its success is subject to the validity of the LK06 algorithm. According to the 
numerical test, the LK06 algorithm suffers approximately 20% errors for the case that the 
number of the member galaxies is not high enough. It implies that the LK06 algorithm is 
definitely restricted to the rich clusters with large numbers of galaxies. Therefore, it may be 
necessary to refine and improve the LK06 algorithm itself for the application to poor cluster 
samples. Our future work is in this direction. 



This work is supported by the research grant No. R01-2005-000-10610-0 from the Basic 
Research Program of the Korea Science and Engineering Foundation. 
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Fig. 1. — Positions of the Virgo member galaxies in the two dimensional projected space. 
The arrow represents the direction of the projected major axes of the Virgo cluster. 
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Fig. 2. — Probability density distributions of the cosines of the angles between the position 
vectors of the Virgo cluster galaxies and four different choices of the Virgo cluster major 
axis in the two dimensional projected plane of sky In each panel, the histogram with 
Poisson errors represents the observational data points from the Virgo Catalog, the solid 
line is the analytic fitting function based on the LK06 reconstruction algorithm, and the 
dotted line corresponds to the case of no correlation. The top left panel corresponds to the 
best-fit result according to which the major axes of the Virgo cluster in the two dimensional 
projected plane is determined, while the other three panels show how the agreements between 
the observational and the analytic results change if different directions other than the major 
axes are used. 
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Fig. 3. — Comparison between the observational and the analytic results for the alignments 
at different two dimensional cut-off radii (i? C ut)- 
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Table 1. The cut-off radius (R cu t) in the two dimensional projected space, the number of 
the member galaxies enclosed within i? cu t, the reconstructed two axial ratios, and the 

best-fit value of the correlation parameter. 



-Rcut 

(Mpc) 


N 9 


a/c 


b/c 


s 


2.24 


1275 


0.53 


0.73 


-0.25 


1.68 


1221 


0.64 


0.69 


-0.25 


1.49 


1154 


0.64 


0.69 


-0.25 


1.12 


902 


0.64 


0.69 


-0.25 



